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ON  THE  INTERPOLATION  OF  GRAVITY  ANOMALIES  AND 
DEFLECTIONS  OF  THE  VERTICAL  IN  MOUNTAINOUS  TERRAIN* 


H.  Baussus  von  Luetzow 
U.S.  Army  Topographic  Laboratories 
Fort  Belvoir,  Virginia  22060 


ABSTRACT :  *  The  paper  first  addresses  the  interpolation  of  gravity  anomalies 
In  mountainous  terrain,  to  be  represented  as  the  sum  of  a  "signal"  variable 
with  a  quasi-stationary  estimation  structure  and  a  computable  "noise"  variable 
without  a  stationary  character.  It  then  develops  the  particular  solution  of  the 
boundary  value  problem  of  physical  geodesy  which  permits  a  similar  representation 
of  deflections  of  the  vertical  and  draws  some  conclusions  concerning  the 
inapplicability  of  Molodensky's  series  approach  and  of  the  collocation  method 
for  the  accurate  determination  of  vertical  deflections  from  unmodified  gravity 
anomalies  in  mountainous  terrain.  Thereafter,  it  discusses  the  estimation  of 
signal -type  deflections  of  the  vertical  by  means  of  spatial  covariance  functions, 
i.e.,  by  a  linear  regression  technique  called  statistical  collocation  in  physical 
geodesy,  and  provides  first  order  expansions  of  planar  covariance  functions. 


1.  INTRODUCTION.  Deflection  of  the  vertical  components  £  and  n  play  a 

* 

role  in  the  adjustment  of  geodetic  networks,  in  the  computation  of  height  anomaly 
differences,  and  in  the  transformation  of  local  coordinates  into  terrestrial 
coordinates.  Short  of  a  three-dimensional  solution  of  the  geodetic  boundary  value 
problem  under  consideration  of  mountainous  terrain,  deflection  components  and  gravity 
anomalies  Ag  are  also  desirable  for  the  numerical  upward  continuation  of  the  first 
order  derivatives  of  the  anomalous  gravity  potential.  The  interpolation  or 
estimation  of  gravity  vector  components  in  flat  terrain  is  not  Inherently  difficult. 
In  mountainous  terrain,  gravity  anomalies  Ag  are  profitably  modified  to  Faye 
anomalies  AgF  by  means  of  terrain  corrections  C,  to  be  followed  by  a  transformation 

♦Paper,  accepted  for  presentation  at  the  Vlth  International  Symposium  on  Geodetic 
Computations,  Munich,  W.  Germany,  Aug  31-Sep  5,  1981.  To  be  published  In  sympo¬ 
sium  proceedings  and  In  Proc.  25th  Conference  of  Army  Mathematicians,  U.S.  Army 
Research  Office,  Research  Triangle  Park,  NC  27709. 
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to  Bouguer  anomalies  Agg  which  permit  an  approximate  two-dimensional  Interpolation. 
Isostatic  gravity  anomalies  Agj,  to  be  corrected  by  the  Indirect  effect,  would 

require  a  three-dimensional  Interpolation  technique  In  the  case  of  high  accuracy.  The 
problem  of  deflection  estimation  has  been  discussed  by  Helskanen  and  Moritz  [1967]  and 
others, Including  the  method  by  Molodensky  et  a*.  [1962J  for  the  calculation  of  deflection 
differences  In  flat  terrain,  and  the  difficulty  to  Interpolate  6,  n  In  rough  moun¬ 
tainous  terrain.  Baussus  von  Luetzow  [1380]  addressed  the  optimal  denslfi cation 
of  deflections  of  the  vertical  In  flat  terrain  with  and  without  consideration  of 
gravity  anomalies  and  extended  Molodensky's  approach.  Badekas  and  Mueller  [1968] 
utilized  Eotvos  torsion  balance  measurement  together  with  appropriate  terrain 
corrections  for  the  Interpolation  of  vertical  deflections,  a  time-consuming  procedure 
and  soon  to  be  replaced  by  the  employment  of  moving  base  gravity  gradiometers . 

Regardless  of  these  efforts,  an  effective  6, n-estlmation  method  applicable  in  mountalnou; 
terrain  will  still  be  valuable  and  may  also  aid  deflection  estimation  under 
consideration  of  a  series  of  discrete  Inertial  measurements.  Section  2  of  this 
study  addresses  the  Interpolation  of  gravity  anomalies  In  mountainous  terrain.  In 
section  3,  the  appropriate  solution  of  the  boundary  value  problem  for  vertical 
deflections  is  presented  and  reformulated  for  optimal  deflection  estimation  of 
"signal"  components  of  £  and  n  and  computation  of  topographic  "noise"  terms.  The 
estimation  of  signal-type  components  by  means  of  spatial  collocation  and  the 
development  of  first  order  approximations  of  spatial  covariance  functions  Is  the 
subject  of  section  4. 

2.  INTERPOLATION  OF  GRAVITY  ANOMALIES.  It  Is  well  known  that  an  accurate 

i 

analytical  representation  of  free-air  anomalies  In  pronounced  mountainous  terrain 
can  only  be  achieved  by  a  polynomial  of  high  degree  by  means  of  Ag-data  available  In 
a  network  of  high  resolution.  As  a  consequence,  satisfactory  linear  Interpolation 
requires  small  mesh  sizes  ax.  Ay.  The  following  modified  anomalies  have  been 
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useful  for  geodetic  applications  and  the  purpose  of  Interpolation: 

Agp  ■  Ag  +  C  (1) 

where  C  Is  the  terrain  correction,  is  called  Faye  anomaly. 

36T 

AgB  ■  gF  -  bh  + -  (2) 

2R 

is  the  modified  Bouguer  or  complete  topographic  anomaly  where  b  «  0.112  mgalnf1 

is  the  Bouguer  gradient,  h  is  the  elevation  of  terrain,  «T  is  the  potential 
of  topographic  masses,  and  R  ■  6371  Km  is  the  earth's  mean  radius. 

A9i  =  A9g  +  Cj  +  a6c  *  Ag  +  C  -  bh  +  Cf  +  a$5  +  r  (3) 

Is  the  Isostatic  anomaly  valid  for  the  compensated  geold  with  a  «  0.3086  mgal  m"1 , 
6z  as  the  vertical  separation  between  geold  and  cogeold,  and  r  as  a  random  error. 
Equation  (3)  nay  be  further  written  as 

Ag  =  Ag,-  +  Ct  +  r  (4) 

where  represents  the  aggregate  of  terms  computable  from  the  known  topography. 

In  a  more  general  formt  also  applicable  to  the  optimal  estimation  of  vertical 
deflections,  equation  (4)  Is  reformulated  as 

m  *  s  +  n  +  r  (5) 


In  this  equation,  m  Is  a  "message"  variable,  s  is  a  "signal"  variable,  n  Is 
deterministic  or  computable  "noise,"  and  r  Is  random-type  noise. 

Under  consideration  of  a  linear  signal  estimation  structure,  a  signal  can 
then  be  optimally  estimated  as 

Se  -  LCot  -  nf  -  r1 )  (6) 


where  L  denotes  a  linear  operator  and  the  subscripts  e  and  1  refer  to  the 


estimation  point  Pe  and  measurement  points  Pf,  respectively.  The  optimal 
measurement  at  Pe  results  as 

me  *  sg  +  ne  +  re  ■  L(ro^  -  nj}  +  re  -  L(r^)  +  ne  (7) 

The  estimation  error  Is 


e(rae)  •  e(se)  +  e[re  -  L(ri)J 


(8) 


The  corresponding  estimation  error  resulting  from  the  utilization  of  topographically 
unmodified  measurements  m^  Is 


e(me)  *  se  -  L(sf)  +  re  -  L(ri)  +  ne  -  L(ni) 
*  e(se)  +  elre  -  L(ri)]  +  ne  -  L(ni ) 


(9) 


I 


Comparison  of  equation  (9)  with  equation  (8)  shows  that  the  non-optimal  Interpolation 
process  Is  associated  with  a  "topographic"  estimation  error  ne  -  L(nf)  which 
becomes  In  general  Intolerable  In  moderate  to  rough  mountainous  terrain  and  thus 
Induces  the  requirement  of  a  fine  mesh  data  grid. 

The  Interpolation  of  Isostatic  anomalies  by  means  of  spatial  collocation  will 
be  treated  In  conjunction  with  the  Interpolation  of  Isostatic  deflections  of  the 
vertical  In  section  4. 

3.  FORMULATION  OF  A  VERTICAL  DEFLECTION  SOLUTION  SUITABLE  FOR  OPTIMAL 
INTERPOLATION.  Gravimetric-topographic  solutions  for  the  anomalous  gravity  potential 
and  deflections  of  the  vertical  which  Inherently  permit  a  "signal-noise"  separation 
according  to  equation  (5)  have  been  established  by  Pelllnen  [1969],  Moritz  [1969]  , 
and  Baussus  von  Luetzow  [1971].  The  latter  emphasized  that  these  essentially 
identical  solutions  are  almost  equivalent  to  those  of  Molodensky  et  al .  [1962]  and  Brovar 
I  19643*  but  are  less  data  dependent,  more  direct  from  the  computational  view,  and 
■ore  advantageous  for  the  utilization  of  artificial  satellite  data.  The  notations  to 
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foe  used  are  the  following: 


t  *  T.y 


r  s  -(-25 — ) 

v3x  'h*const. 

n  “  '^Sy- \i-const. 
h  =  h* 


81 B arc  It 

v  ~  It 

d  a 

'  3y 


k  «  6.67.10“8cm3g-1sec"2 
p  *  2.67  g  cm-3 
R  »  6371  Km 
10  *  2R  sin  | 

1  *  (rA  +  rp  ‘  2rArP  C0S  ^ 


Ag  »  g-Y 


A9p  ■  Ag  ♦  C 


height  anomaly 

anomalous  gravity  potential 

normal  gravity 

prime  vertical  deflection 

meridian  vertical  deflection 

elevation  of  terrain  referring  to  moving  point  P 

elevation  of  terrain  referring  to  fixed  computation 
point  A 

northern  terrain  inclination 

eastern  terrain  inclination 

derivatives  taken  along  the  local  horizon  in  a 
northern  and  eastern  direction 

global  mean  gravity 

azimuth  angle  counted  clockwise  from  north 

angle  between  the  radius  vectors  rA  and  rp 
originating  at  the  earth's  spherical  center 

Stokes'  function 

gravitational  constant 

standard  density 

earth's  mean  radius 

see  Figure  1 

see  Figure  1 

unit  sphere  (full  solid  angle) 
measured  gravity 
gravity  anomaly 
terrain  correction 
Faye  anomaly 
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Bouguer  gradient 

potential  of  topographic  masses 


b  *  0.112  mgal  m-1 
6T 

ag«  *  a9f  "  bh  +  modified  Bouguer  or  complete  topographic 

*K  anomaly 

The  geometry  Involving  hA,  h=hp,  ♦,  R,  10,  and  1  Is  evident  from  Figure  1  below. 

A 


Figure  1 


The  established  first  order  solution  for  the  deflection  components  Is 
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It  should  be  noted  that  s g  is  in  general  very  small  and  that  the  computation 
of  6g1(6c»  and  6n  requires  only  Integrations  over  oi«  4».  It  has  further  to 
be  emphasized  that,  according  to  Baarda  [1979],  the  Inclination  angles  and 
$2  should  not  exceed  70°. 

Equation  (10)  Is  now  reformulated  under  consideration  of 

Agp  *  Ag  +  (Agp  -  Ag)  »  Ag  ♦  6g2  (14) 

AgB  s  Ag  +  (AgB  -  Ag)  *  Ag  +  6g3  (15) 

In  these  equations,  Ag  Is  a  signal  variable,  profitably  the  Isostatic  anomaly 
defined  In  equation  (3).  In  comparison  with  6g2,  6g3  Is  a  relatively  smooth 

topographic  quantity. 

The  substitutions  (14)  and  (15)  transform  equation  (10)  into 

A  R  CCrKn  +  r  f A nM  icos“l  ds(^-)  da  -  Ag  +  Gi(Ag)  Itn  Bi) 
hrttJJW  +  Gi  (Ag)J  isinaf  dtp  °°  G  (tn  e2f 

a 

+&/f\ [«fc  *  «* ♦ «.(•«,»]  (sa^* ♦ a  (sa  061 

a 

The  first  two  terms  of  equation  (16),  Involving  the  anomaly  Ag,  represent  the 
"signal"  components  of  K  and  n.  The  following  three  terms  constitute  computable 
topographic  "noise."  Permitting  for  random-type  errors  r^  and  r^,  equation  (16) 
can  be  written  in  analogy  with  equations  (4)  and  (5)  as 

{SI  *  (,7) 

The  numerical  determination  of  the  three  topographic  terms  of  equation  (16) 

Is  a  complex  task,  which  can,  however,  be  accomplished  without  Inherent 
difficulties  by  means  of  high-speed  computers.  In  this  respect,  the  Integration 
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area  relating  to  the  first  topographic  term  can  be  considerably  restricted, 
appears  that  the  last  two  topographic  terms  are  particularly  subject  to  rapid 
changes  in  mountainous  terrain.  Accurate  Interpolation  Is  further  favored  if 
given  and  estimated  deflections  refer  to  points  associated  with  small  terrain 
Inclinations. 

In  accordance  with  Moritz  [1969  ],  the  second  order  correction  for  the 
height  anomaly  is 

6;(2)  s  ff  +  693)  s^da  "  ff  (A9  +  $93)  da 


where 


G2  =  Ma9  +  693  )da  +  (A9  +  693)  tn2Bm 


Here,  Bm  represents  the  maximal  terrain  inclination. 
The  second  order  deflection  corrections  are  then 


J6e(2H  _  (6c(2)(A9))  R  ffr  /  •  ,  j  COSa  1  dS(<|/) 

"  W2>(*9 )f  +  |  sina )  d*  da 


+  3Ri  f A*  (h-hAl2  1  cosal  £jj.  da 
4ir  JJ693  10“  1  sina  I  d* 


Designatlng  the  Integral  terms  of  equation  (20)  as  second  order  topographic 
corrections  and  6nf2\  equation  (17)  assumes  the  modified  form 


o-dsi  •  ci  •  (Si*  csi 


Higher  order  topographic  correction  terms  are  not  warranted  because 
of  a  decreasing  convergence  radius  In  connection  with  higher  derivatives,  the 
assumption  of  a  standard  density  or  density  uncertainties,  respectively,  and 


imperfect  isostatic  equilibrium.  The  structure  of  equation  (21)  clearly  exhibits 

the  fact  that  a  highly  accurate  computation  of  and  t/2)  cannot  be  achieved 

by  the  exclusive  utilization  of  free  air  anomalies  Ag.  For  the  same  reason,  iterative 

solutions  of  the  integral  equations  for  generalized  surface  densities  by  Molodensky 

et  al .  [1962]  and  Brovar  [1964]  and  the  series  solution  by  Molodensky  et  al .  [1962] 

in  general  do  not  converge  in  mountainous  terrain.  The  latter  permits  for 
auxiliary  boundary  surfaces  under  utilization  of  a  shrinking  parameter  k  >*0  and  thus 

implies  the  possibility  of  analytical  continuation  with  p=0.  For  the  same  reason, 
collocation  solutions  would  only  satisfactorily  apply  with  respect  to  signal  variables 
l  and  fi.  The  analytical  upward  continuation  jof  the. first  derivatives  of  the  anomalous 
gravity  potential  in  mountainous  terrain  would  require  a  supplemental  approach. 

4.  SIGNAL  ESTIMATION  BY  STATISTICAL  COLLOCATION  AND  FIRST  ORDER  EXPANSIONS 
OF  PLANAR  COVARIANCE  FUNCTIONS.  As  indicated  by  Baussus  von  Luetzow  [1980], 
deflection  differences  in  flat  terrain  may  be  advantageously  determined  by  a 
combination  of  statistical  collocation  and  Vening  Meinesz  formulae  provided 
gravity  anomalies  are  also  available  in  sufficient  density  within  a  limited 
region.  Four  point  deflection  estimation  errors  with  mesh  sizes  Ax  =  5  km, 

B  km,  and  24  km  were  found  to  be,  respectively,  of  the  order  0.1  arcsec,  0.2 
arcsec,  and  1.0  arcsec  in  the  case  of  estimators  free  of  errors.  Astrogeodetically 
determined  deflections  are,  however,  presently  associated  with  errors  of  the 
order  of  0.25  arcsec.  In  accordance  herewith,  it  is  advantageous  to  employ  a 
relatively  great  number  of  estimators  if  this  is  feasible. 

A  0^ 

The  signal  variable  to  be  estimated  and  representing  either  £  or  n  may  be 
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x  ,  and  the  estimators  may  be  written  x±  with  6A  as  a  correlated  measurement 

* 

error  independent  of  Under  the  assumption  of  an  existing  signal  and  noise 
covariance  structure  the  following  linear  regression  equation  can  be  formulated: 

xe  +  6i)  *  +Ai)  (22) 

It  is  then  in  matrix  form,  with  bars  indicating  covariances. 


Ai(xii 


aA> 


Vik 


1.2. 


(23) 


The  solution  for  the  regression  coefficient  matrix  follows  as 

Ai  =  *e  N1k  1  (24) 

In  the  case  of  given  astrogeodetic  vertical  deflections,  fij  may  be  composed 

2 

of  astrogeodetic  errors  with  a  variance  (0.25  arc sec)  and  a  correlated  error 

partially  caused  by  imperfect  isostatic  equilibrium. 

VJith  respect  to  the  basis  for  the  statistical  collocation  approach  in 

physical  geodesy,  reference  is  made  to  Bjerhammar  [1973],  Grafarend  [1973], 

Krarup  [1969],  Lauritzen  [1971],  Moritz  [1970],  and  Tscherning  [1973].  Of 
significance  is  that  the  spatial  covariance  function  for  the  disturbing  gravity 
potential  has  to  satisfy  Laplace's  equation.  Baussus  von  Luetzow  [1973] 
emphasized  the  necessity  to  treat  c-if  as  a  correlated  random  variable  where  c  is 
a  deterministic  development  of  c  in  spherical  harmonics  of  at  least  degree  and 
order  15.  In  accordance  herewith,  the  requirement  of  homogeneity  prescribes  and 
at  least  permits  in  practice  a  restriction  to  the  planar  approach  in  physical 
geodesy.  Accordingly,  |j  *  -Ag ,  |^,  and  |y  are  supposed  to  satisfy  Laplace's 
equation.  It  Is  realized  that  the  convenient  requirements  of  homogeneity 
and  quasi-flat  terrain  are  only  approximately  satisfied. 
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Moritz  [1976]  and  Nash  and  Jordan  [1978]  established  specific  T-covarlance 
functions  which  can  be  expanded  Into  space  In  a  closed  form.  As  has  been 
shown  by  the  latter  authors,  the  spatial  covariance  function  Is 


♦TT(r,  Zi,  22)  =  y*  pHp)®  1  2  Oo(rp)dP 

O 

00 

F(P)  *  y*  r$TT(r)J0(rp)dr 

O 


(25) 

(26) 


In  equations  (25)  and  (26),  j0  Is  the  zero-order  Bessel  function,  r  Is  the 
variable  planar  distance,  zj  and  Z2  are  the  elevations  of  two  points,  and 

^tt  (r)  t^ie  P^anar  T  "  covariance  function. 

The  spatial  vertical  deflection  covariances  may  be  derived  from 
equation  (25)  in  the  form 


(*a  (r.  zi,z2)  1 

=  *  (w-rfl  ®TT^r*  zi»  z2 ) 

Unr|  (r,  Zi  ,z2) " 

Sy2 

where  Vi  sY(zi),ir2  »y(z2). 

Forij^-  functions  which  permit  the  derivation  of  realistic  vertical  deflection 
covariance  functions,  the  Hankel  transforms  (26)  and  (25)  cannot  be  evaluated 
in  closed  form.  As  an  example,  Jordan’s  [1972]  third-order  Markov  model 


(27) 


(28) 


leads  to  the  hypergeometric  function  when  introduced  in  equation  (26).  Thereafter, 
♦^(r,  Zi ,  z.2)  only  can  be  obtained  by  an  extremely  lengthy  numerical  integration. 
For  this  reason.  It  appears  to  be  advantageous  to  develop  first  order  approximations 
for  spatial  vortical  deflection  covariance  functions  under  consideration  of  Jordan's 
[1972]  planar  results.  In  this  respect  It  has  to  be  emphasized  that  Jordan 
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Interchanged  the  conventional  partial  differentiations 


Figure  2 

Under  consideration  of  Figure  2,  It  Is  to  the  first  order 

«1  *  «0  ♦  i  h  ‘  E«  * 

It  is  then,  to  the  first  order. 


*  V,  HC  +  5.1“*,  *  wgw*,  ♦  V 


lis.7 


Mf 


It  is  further  at  level  2  =  0 


\  si  MS. 

3z  ”  ”  3y  3z  ”  ”  ay  '  G  7  G  ay 
so  that 

r  He  -  1  3Age  l  a_— — 
to  32  E  ay  £  3y  ^oA9e 

Under  consideration  of 

toAge  ‘  "  5?»o  =  “  h<r>  f 

It  Is 


(29) 


(30) 


(31) 


(32) 


(33) 


|y  [h(r)f  ]  -  ^+h(r)(~:  -  £r)  *  §7  cos2o  +  7  s1n2a 


(34) 


The  final  results  are  hereafter 
__  3oran  h  3h 

5^2  (»*,  O,  Zj,  z2)  ”  vf^G  ^  s1n2<x  +  9r  cos2a^zi  +  Z2*  ^ 

3o  oh  l 

^  (r,  a,  z1$  z2)s*nn-^|(r  cos2°+77s1n2a^zi  +  z 2)  (36) 

where  ^  and  $nn  represent  the  planar  covariance  functions  and  nDne } 

respectively,  and  where  Og«  o^g-cvarag)15. 

In  analogy  with  equation  (30),  It  is 

-9i  4g2  =  Ag0  age  +  ag0  ||9e  (Zi  +  z2)  (37) 

Under  utilization  of  the  planar  approximation 


3Aq 

3Z 


.6(iL  +  lH) 
v3y  3x' 


(38) 


it  is 


-GA9n( 


3  He 


‘3y  3x  ^  ”^3y  gx  A9one) 

With  the  aid  of  equation  (33),  equation  (39)  can  be  formulated  as 


(39) 


3£e  3a  oa  h  3h 

-GAgft(—  +  —  )=  -G  -£-2-  (-  +  — ) 
y0v3y  ax  '  y/2  'r  3r' 


(40) 


Accordingly, 


49.  4gw  (r,  z, 


z2^  “  *gg 


1^1  (**  +  »!*)  (2  +  z  ) 
y/f-  'r  3  r'  1  l  2 7 


(41) 


where  *gg  -  Ag0*ge. 

It  is  evident  from  equations  (35),  (36),  and  (41)  that  these  represent  convenient 
closed  approximations  of  the  three  spatial  covariance  functions  of  particular 
Interest.  In  general,  the  planar  covariance  functions  should  apply  to  the  lowest 
z-level  In  a  particular  area  of  application. 
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5.  CONCLUSION.  The  Immediate  accurate  interpolation  of  gravity  anomalies 
and  deflections  of  the  vertical  In  mountainous  terrain  is  only  possible  from 
data  provided  In  a  grid  of  high  resolution.  Optimal  Interpolation  from  data 
given  at  points  separated  by  distances  of  the  order  5-10  km  or  from  multiple 
data  Incorporating  measurement  noise  with  shorter  spacing  can  be  accomplished 
by  an  appropriate  representation  of  gravity  anomalies  and  deflections  as  a 
signal -noise  process  with  nonstationary  noise  computable  from  the  earth's 
topography.  In  the  case  of  deflections,  a  special  solution  of  the  geodetic 
boundary  value  problem  Is  required.  As  a  first  approximation,  Faye  anomalies 
may  be  used  as  signal  variables.  Isostatic  anomalies  modified  by  the  indirect 
effect  provide  a  greater  degree  of  homogeneity  and  isotropy.  Implementation 
of  the  theory  requires  the  utilization  of  existing  "isostatic"  computer  programs 
and  the  establishment  of  a  supplemental  program  under  consideration  of 
furnished  analytical  solutions.  Signal  estimation  has  to  be  facilitated  by 
the  use  of  spatial  covariance  functions  first  order  approximations  of  which 
may  be  computed  relatively  easy  from  planar  covariance  expressions.  The 
optimal  interpolation  method  in  conjunction  with  the  special  solution  for 
deflections  Indicates  that  Iterative  or  series  solutions  of  the  boundary 
value  problem  of  physical  geodesy  cannot  be  expected  to  converge  in  mountainous 
terrain.  The  method  developed  Is  of  practical  significance  for  the 
densificatlon  of  gravity  anomaly  and  astrogeodetic  deflection  networks  in 
mountainous  terrain  and  is  also  valuable  or  Indispensable,  respectively,  for 
the  optimal  estimation  of  gravity  anomalies  and  deflections  from  astrogeodetic 
and  Inertial  data  in  mountainous  areas. 
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